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First AU-scale observations of V1647 Ori witli VLTI/IVIIDI 
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Abstract. The young eruptive star VI 647 Ori was observed with MIDI, the mid-infrared interferometric instrument at the 
Very Large Telescope Interferometer (VLTI), on March 2, 2005. We present the first spectrally resolved interferometric visi- 
bility points for this object. Our results show that (1) the mid-infrared emitting region is extended, having a size of ~7 AU at 
10 /im; (2) no signatures of a close companion can be seen; (3) the 8 — 13 ^m spectrum exhibits no obvious spectral features. 
Comparison with similar observations of Herbig Ae stars suggests that VI 647 Ori probably possesses a disk of moderate flar- 
ing. A simple disk model with T^r~^'^"^, E^r"^'^, Ma ~ 0.05 Mq is able to fit both the spectral energy distribution and the 
observed visibility values simultaneously. 
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1. Introduction 



techniques: interferometric - stars: circumstellar matter - stars: individual: V1647 Ori - in- 



In January 2004 a new reflection nebula (McNeil's Nebula) 
appeared in the LDN 1640 dark cloud of the Orion B molec- 
ular cloud complex (" McNeil et al. 20 04t. V16470ri, whose 
outburst (~4mag in the I-band) caused the appearance of 
McNeil's Nebula, is a low-mass pre-main sequence object 
JBriceno et al. 20041 [Abraham et al. 2004a' . Its eruptive be- 
haviour suggests that VI 647 Ori is either an FU Orionis (FUor) 
or an EX Lupi (EXor) type object, or maybe an intermediate- 
type object between FUors and EXors (Muzer oUe et al. 20051 
[Kospal et al. 2005 1. Near-infrared colour maps show that the 
source is embedded in an elongated disk-like structure, whose 
size is approximately 7000 AU, and its inclination (the angle 
between the normal of the disk and the line of sight) is about 
60° (Acosta-Pulido et al., in prep.). The object had been gradu- 
ally fading until October 2005, when the eruption rapidly ended 
HKospal et al.^OOS) . 

Both FUors and EXors are low-mass pre-main sequence 
stars which exhibit optical brightening of several mag- 
nitudes. FUors are characterised by outbursts of several 
decades, while EXors exhibit repetitive outbursts on monthly 
timescale. The eruptive mechanism of both types is thought 
to be a rapid temporal increase of the disk accretion rate 
l |Hartmann & Kenyon 1996| and references therein). In this 
model the young star is accreting a substantial amount of ma- 
terial from the parent molecular cloud core (10^^ Mq/yi) via 
a circumstellar disk. The infalling matter piles up in the in- 
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ner disk until its surface density - and thus opacity - be- 
comes high enough to switch on a thermal instability lead- 
ing to the dramatically increased accretion rate (10^^ - 10~^ 
M0/yr). It is still debated whether perturbation due to a close 
companion is needed to trigger the outburst, and it is also an 
open question whether all eruptive stars have companions (e.g. 
IWang et al. 2004||keipurth & Aspin 2004b> . 

The geometrical structure of the inner part of the cir- 
cumstellar material in eruptive systems - though it is prob- 
ably closely related to the outburst mechanism - is still un- 
clear. Most models assume that apart from the circumstellar 
disk, there is an infalling remnant envelope of the molecular 
cloud core, whose inner radius is of the order of a few AU 
( |Kenyon & Hartmann 1 99 1| [Turner et al. 1997> . This envelope 
supplies the disk with the material needed for the strong erup- 
tions. In the case of V16470ri, Muzerolle et al. (2005) pro- 
posed that an optically thin envelope would be necessary to 
explain the flat infrared spectral energy distribution. However, 
a direct detection of such envelopes around eruptive stars has 
not been attempted so far. 



In this paper we report on AU-scale observations of 
VI 647 Ori at 10 fim with MIDI, the mid-infrared interfer- 
ometric instrument mounted at the Very Large Telescope 
Interferometer (VLTI) of ESO's Paranal Observatory 
(ILeinert et al. 20031 . We present spectrally resolved inter- 
ferometric visibility points, and analyse the data with special 
emphasis on the structure of the circumstellar material and on 
the signature of a nearby companion. 



2 



P. Abraham et al.: VLTI/MIDI observations of VI 647 Ori 



2. Observations and data reduction 

V1647 Ori was succesfuUy observed with MIDI on the UT3- 
UT4 baseUne of the VLTI on March 2, 2005. The projected 
baseline length was 56m with a PA=112°. Due to the lack of 
adequate guide star, MACAO', could not support the obser- 
vation, which was thus affected by the seeing (0.9-1.0" at the 
time of the observation). Single telescope N-band spectra were 
also taken with both UTS and UT4. Additional N-band spectra 
were obtained on December 3 1 , 2004 (then the interferometric 
part of the observation failed). On March 2, HD 37160 was ob- 
served as a calibrator, while on December 31 observations of 
HD 50778, performed 2 hours later, were used. 

The obtained data set consists of acquisition images with 
the N8.7 filter, 8-13 /im low resolution spectra (R=30), and 
interferometric measurements. In the data reduction we fol- 
lowed the general processing scheme as described by Leinert 
et al. (2004). The measured raw visibilities of V1647 Ori were 
divided by raw visibilities obtained from the HD 37160 data. 
The observations were reduced in two independent ways: with 
the MIDI Interactive Analysis (MIA) package, which uses 
the power spectrum method, and the Expert Work Station 
(EWS) package, which is based on a coherent, linear averag- 
ing method. For the extraction of the spectra the two packages^ 
consider fitted and fix masks, respectively. The interferometric 
results were then compared and showed agreement better than 
10%. 

All available calibrator data obtained on March 2, 2005 
were reduced in order to determine the errors of the spectrally 
resolved visibilities. After investigating the data we decided to 
apply a uniform multiplicative error of 10% as a conservative 
estimate. 

At both epochs the N-band spectra obtained by the two tele- 
scopes differ in their absolute flux levels by a factor of « 1 .4, 
and their spectral shapes slightly deviate at A>12 /im. The rea- 
son of this artifact is not yet clarified, thus we evaluated the 
acquisition images and extracted flux density values at 8.7 /im 
for VI 647 Ori (for the calibrator, HD 37160, 1 1.85 Jy was as- 
sumed). Photometry taken with UT3 turned out to be more re- 
liable, judged from the repeatibility of the standard star mea- 
surements, and we estimate its absolute accuracy to be 15%. 
The absolute calibration was then carried out by scaling the N- 
band spectra to the photometric points at 8.7 /tm. 

3. Results 

Figure[l]shows the calibrated visibilities as a function of wave- 
length. The values gradually decrease from 0.87 to 0.75 be- 
tween 8 and 10 /xm, and are approximately constant at longer 
wavelengths, with a small dip between 12 and 13 /im. This pat- 
tern clearly demonstrates that the emitting region is resolved, 
and suggests a non-uniform temperature distribution of the 
emitting material (otherwise the visibilities would monotoni- 
cally increase towards longer wavelengths due to the lowering 
spatial resolution). 

' http://www.eso.org/projects/aot/macao_vlti/, the Curvature AO 
system for the VLTI 
^ http://www.mpia-hd.mpg.de/MIDISOFT/ 
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Fig. 1. Calibrated visibilities as a function of wavelength. The 
uniform error bars of 10% reflect our conservative estimate of 
the uncertainties. The solid line represents the model visibility 
curve of V1647 Ori (see Sect 4.). 



Adopting Eq. 1 from Leinert et al. (2004), we calculated 
visibility values for a series of extended model sources of 
increasing sizes, represented by Gaussian brightness profiles, 
and compared them with our observations at different wave- 
lengths. The derived FWHMs increase towards longer wave- 
lengths from 6^3 mas at 8.2 /im to 16±3 mas at 13 fim. At the 
distance of VI 647 Ori (450 pc) these angular sizes approxi- 
mately correspond to 2.7 and 7.2 AU linear sizes, respectively. 
The increasing sizes of Gaussians with increasing wavelength 
indicate a radially decreasing temperature profile. 

Fig. 12 presents our 8-13 //m spectrum of V1647 0ri ob- 
tained with MIDI on UT3 in March 2005. The spectrum is 
smoothly rising towards longer wavelengths, and no obvious 
spectral features (silicates, PAHs, molecular ices) can be seen. 
The shape of the December 2004 MIDI spectrum is identical to 
the plotted spectrum within the measurement uncertainties. For 
comparison, in Fig.|2]we also overplotted an earlier spectrum of 
V1647 0ri from March 2004 (see figure caption). Comparison 
of the spectra suggests that the N-band spectrum changed dra- 
matically (especially at longer wavelengths) between March 
and December 2004, but exhibited an approximately constant 
spectral shape afterwards, between December 2004 and March 
2005. 

We calculated the spectrum of the correlated flux as the 
product of the measured N-band spectrum and the spectrally 
resolved visibilities, in order to study the emission of the inner- 
most part of the circumstellar structure. We found that the main 
fraction of the observed N-band flux (w 70%) is emitted in the 
inner regions. The relatively flat shape of the correlated spec- 
trum indicates that the innermost part has higher temperature. 
Like in the full spectrum, no spectral features can be observed 
in the innermost spectrum. 

4. Discussion 

V1647 0ri is a l ow-mass object JBriceiio et al. 20041 
lAbraham et al. 2004ai . Nevertheless, due to its increased 
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Fig. 2. N-band spectra of V1647 Ori at different epochs: solid 
line - VLTI/MIDI (March 2, 2005, this study); dotted line 
- UKIRT/Michelle (March 11, 2004, Andrews et al. 2004). 
The MIDI spectrum was scaled to the photometric point at 
8.7 ^m obtained with UT3 {filled dot). Neither the MIDI nor 
the UKIRT spectrum shows any obvious spectral features. 



outburst luminosity (Lboi — 4:ALq, IMuzeroUe et al. 2005> . 
it may look similar to intermediate-mass young stars. As 
the first step of our analysis, we compare our MIDI results 
with similar observations of a sample of Herbig Ae/Be stars 
JLeinert et al. 2004 1. Looking first at the visibility curves 
(Fig-E of the present paper; Fig. 4 in Leinert et al. 2004) 
one notices a general agreement between the shapes of the 
V1647 0ri and the Herbig Ae/Be curves: a faster drop in 
the 8 — 9 fim range is followed by almost constant visibility 
values at longer wavelengths (exceptions are HD 179218 
and 510ph). All Herbig Ae/Be stars show lower visibility 
values than VI 647 Ori, indicating that these objects are more 
resolved. This finding can be explained to a large extent by 
distance differences: these Herbig Ae/Be stars are at a distance 
of 103 - 244 pc compared with the distance of VI 647 Ori 
of 450 pc. We tested this hypothesis with the help of a disk 
model (see below) in which the distance was scaled from 450 
to 150pc, and found visibility curves similar to those of the 
Herbig Ae/Be stars. 

In Fig.|3l we reproduced Fig. 5 of Leinert et al. (2004). 
In order to overplot VI 647 Ori in the diagram, we com- 
puted the half-light radius as defined in Leinert et al. (2004), 
and estimated the [12] — [25] magnitude difference by inter- 
polating among photometric points measured by the Spitzer 
Space Telescope in March 2004 (IMuzerolle et al. 200 5». When 
plotting the Herbig Ae/Be stars, different symbols were 
used for Group I sources (thought to harbour flared disks, 
IMeeus et al. 20011 . and Group II sources (believed to possess 
non-flaring disks). VI 647 Ori fits into the trend that redder 
disks are more extended. Considering its position between 
Group I and Group II sources one may conclude that V 1 647 Ori 
probably has a disk of moderate flaring. However, more com- 
plex models including additional components, such as en- 
velopes, cannot be excluded. 
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Fig. 3. CoiTelation between the mid-infrared spectral shape, 
calculated as [12] - [25] = -2.5 log(i^j.(12 ^m)/F,,(25 /im)), 
and the half-light radius (adapted from Leinert et al. 2004, see 
text). V1647 0ri is marked with a filled circle. 



FU Orionis, the prototype of the FUor class, was also ob- 
served with MIDI ( Quanz et al. 2006 1, and - being located 
at the same distance - it is natural to compare FU Ori and 
VI 647 Ori. The spectrally resolved visibility curves of FUOri 
are very similar to that of V 1 647 Ori (at least for 2 of the 3 mea- 
sured baselines) in both wavelength dependence and absolute 
value. Their N-band spectra are, however, somewhat different: 
FUOri exhibit weak silicate emission, while VI 647 Ori shows 
only smooth continuum radiation. Malbet et al. (2005) mod- 
elled FUOri with a geometrically thin accretion disk which 
could reproduce their near-infrared interferometric results as 
well as the broad infrared-millimetre SED. The same model, 
however, failed to match the mid-infrared visibility curve ob- 
tained by MIDI ( IQuanz et al. 2006> . 

In order to model the spectrally resolved visibilities of 
VI 647 Ori, first we constructed a spectral energy distribu- 
tion. Fortunately, many observations were performed by dif- 
ferent groups in February-March 2004, which cover the whole 
0.5 — 850 /im range and are now available in the literature. 
Though this date was one year earlier than our MIDI observa- 
tion, VI 647 Ori was in the outburst phase at both epochs. The 
SED is plotted in Fig. |3 (for references of individual observa- 
tions see figure caption). The SED is characterised by a flat 
spectral shape (t/_FiyWconst.) in the 3 — 30 /im range, similarly 
to a number of FU Orionis objects, like V1057Cyg, V346Nor, 
ZCMa ( Abraham etal. 2004b, Miflan-Gabet et al. 2005 1. On 
the other hand, the SED of VI 647 Ori significantly differs from 
that of FU Ori dMalbet et al. 2005i . 

We model the SED using a simple disk model. More 
complex configurations, e.g. flat disk plus envelope, might 
also be plausible, but the lack of any silicate emission in 
the N-band spectrum (Fig. |2} argues against the presence 
of an optically thin envelope (e.g. like the one proposed by 
IMuzerolle et al. 20^ ). Our analysis tests the adequacy of sim- 
ple disk models to simultaneously fit the SED and the MIDI 
visibilities. In our model one can specify the radial profiles of 
the temperature and the surface density in the disk, assuming 
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that both follow a single power-law. We assumed surface den- 
sity as E'^r^^'^, an extinction of Ay — 10 mag, and an incli- 
nation angle of 60° (Acosta-Pulido et al., in prep.). The central 
object was an FOV star, but its detailed properties did not af- 
fect considerably the fit. The fitted model parameters for the 
disk were the following: temperature T{1 AU) — 680 K and 
T~r~°'^^, inner and outer disk radii 7Rq and 100 AU, respec- 
tively, disk mass Aid — 0.05 Mq. The best fit to the spectral en- 
ergy distribution is overplotted in Fig.|4] The simple disk model 
matches all data from 0.5 to 3000 /im with high precision. 

We predicted visibilities from the above disk model ac- 
cording to the van Cittert-Zernicke theorem, i.e. the complex 
visibility is computed as the Fourier transform of the source 
brightness distribution, and vice versa. The resulting spectrally 
resolved visibility curve for the exact projected baseline of 
the V1647 0ri observation is overplotted in Fig.^ The pre- 
diction agrees with the observed visibility points within the 
measurement uncertainties. Thus our simple disk model of 
VI 647 Ori fits the spectral energy distribution and matches 
the visibility curve simultaneously. The model is also able to 
reproduce the slope of our N-band spectrum. We note, that 
the canonical model usually proposed for FUors (an opti- 
cally thick, geometrically thin accretion disk with T^r^''-^^ 
(Shakura & Sunyaev 1973), or a flared disk with T^r"*^^^ in 
the inner disk and a shallower temperature profile in the outer 
regions) does not fit the data. However, our model with a shal- 
lower temperature profile of T^r^"^^, whose exponent is sim- 
ilar to that of flared disk regions in standard models, can suc- 
cessfully reproduce the observations. Also, based on recent 
near-infrared interferometric measurements, Millan-Gabet et 
al. (2005) found that many FUors require a different model than 
the canonical one. An exception is FU Orionis itself, which was 
successfully modelled by Malbet et al. (2005) using a geomet- 
rically thin accretion disk. 

Temporal evolution of the N-band spectrum of VI 647 Ori 
(Fig.|2} may suggest a more complex picture of the circumstel- 
lar structure. The different fading rates at 8 and 13 /im might 
indicate the presence of two physically separate components, 
which have different energy budget and cooling rate. Detailed 
time-dependent models will be necessary to describe simulta- 
neouly the SED, the visibilities and the temporal evolution. 

A companion would cause sinusoidal variations in the spec- 
trally resolved visibilities with appropriate baseline position 
angles. The shape of our visibility curve suggests that no com- 
panion is present at the measured position angle whose sepa- 
ration is less than 100 AU and brightness ratio is greater than 
10%. Nor do the acquisition images show any companion. The 
presence of companions is an important issue to discriminate 
between different hypotheses explaining the outburst mecha- 
nism ( |Reipurth & Aspin 2004 b ). Due to the lack of more MIDI 
observations on baselines of different length and orientation, 
we can exclude the companion only along the measured base- 
line. The recent rapid fading of VI 647 Ori ( Ko spal et al. 2005j l 
may prevent us from obtaining more visibility data in the 
present outburst period. A future brightening might offer the 
possibility to clarify the existence of a companion of bright- 
ness somewhat similar to the primary. 
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Fig. 4. Spectral energy distribution of V1647 Ori. All opti- 
cal and infrared data plotted were obtained in February-March 
2004. The thick line is the best fitting disk model (see text); 
dotted line represents contribution from the FOV star Sources 
of data: VRI and JHK: Acosta-Pulido et al., in prep.; g'r'i' 
and JHK: Reipurth & Aspin 2004a; 3.6-70 fim: MuzeroUe 
et al. 2005; 3.8, 4.8 jim: Vacca et al. 2004; submm and mm- 
wavelength: Andrews et al. (2004), Lis et al. (1999), Mitchell 
et al. (2001), Tsukagoshi et al. 2005. 
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